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DIFFERENTIAL  ABSORPTION  LIDAR:  EFFECTS  OF 

SPECKLE  NOISE 

KJELL  OSTBERG 


1.  INTRODUCTION 

There  is  today  a  need  for  methods  for  remote  monitoring 
in  real  time  of  air  pollutants.  Among  the  methods  studied 
for  this  purpose  during  the  last  years,  the  differential 
absorption  lidar  (DIAL)  technique  has  been  found  to  be 
the  most  promising  one  ^(Kildal  &  Byer,  1971;  Byer  & 
Garbuny,  1973;  Hinkley,  1976). 

cIn  the  DIAL  technique  one  uses  a  laser  radar  with  a 
frequency  which  can  be  tuned  in  the  neighborhood  of  an 
absorption  line  of  a  gas  in  the  atmosphere.  One  measures 
the  intensity  of  light  backscattered  from  atmospheric 
aerosol  particles,  a  topographical  target,  or  a  retrore- 
flector.  From  the  ratio  of  backscattered  intensity  at  two 
wavelengths,  one  on  and  one  off  the  particular  absorption 
line,  the  gas  content  in  the  air  column  between  the  laser 
radar  and  the  scattering  volume  can  be  obtained.  * —  | 

It  is  advantageous  to  make  DIAL  measurements  in  the 
infrared  region  of  the  spectrum.  Most  molecules  have 
absorption  lines  in  that  region,  and,  furthermore,  infrared 
lasers  are  eye-safe.  However,  the  direct-detection  tech¬ 
nique,  which  has  mostly  been  used  for  DIAL  measure¬ 
ments  so  far,  is  comparatively  insensitive  in  the  infrared. 
The  sensitivity  of  the  DIAL  measurement  will  be  limited 
by  the  thermal  noise  in  the  detector.  It  was  therefore 
suggested  by  Inaba  &  Kobayasi  (1976)  and  Kobayasi  & 
Inaba  (1976)  that  one  should  instead  use  heterodyne  de¬ 
tection  in  the  DIAL  scheme.  Then  it  would  ideally  be 
possible  to  obtain  quantum -limited  operation,  i.e.,  shot 
noise  induced  by  the  local  oscillator  would  become  the 
dominating  noise  source.  They  estimated  that  the  DIAL 
sensitivity  could  be  increased  by  several  orders  of  magni¬ 
tude.  Later,  Menzies  A  Shumate  (1976)  and  Menzies 
(1978)  made  successful  DIAL  measurements  with  hetero¬ 
dyne  detection  of  ozone,  nitric  oxide,  and  ethylene.  In 
these  measurements  they  used  CW  lasers,  and  the  light 
was  backscattered  from  a  retroreflector  or  from  rough 
surfaces. 

With  a  pulsed  lidar  system  and  utilizing  the  light  back- 
scattered  from  atmospheric  aerosol  particles,  it  is  possible 
to  obtain  range-resolved  measurements.  However,  if  het¬ 
erodyne  detection  is  used  in  this  case,  there  will  be  large 


fluctuations  from  pulse  to  pulse  in  the  received  signal. 
The  reason  is  that  the  backscattered  electromagnetic 
fields  from  all  the  aerosol  particles  interfere  with  each 
other  to  set  up  the  total  backscattered  field.  In  a  time  of 
the  order  of  A/r<*  1 ps  (A — wavelength,  v  =  typical  aerosol- 
particle  speed),  the  particles  have  reshuffled  so  that  a  new 
sample  of  the  random  backscattered  field  is  obtained.  If 
the  pulse  time  tp  is  <A/v ,  the  coherence  time  instead  is 
set  by  since  after  that  time,  there  also  is  a  new  inde¬ 
pendent  sample  of  the  backscattered  field.  This  phenom¬ 
enon  has  been  studied  for  a  long  time  in  connection  with 
microwave  radars  (see,  e.g.,  Marshall  &  Hitschfeld,  1963; 
Wallace,  1963).  In  optics  these  fluctuations  in  received 
intensity  are  named  speckle  noise  (Dainty,  1975;  Good¬ 
man,  1976).  It  is  to  be  expected  that  this  speckle  noise 
will  seriously  degrade  the  sensitivity  of  a  range- resolving 
DIAL  system  with  heterodyne  detection.  The  operation 
of  the  system  will  then  no  longer  be  quantum  limited. 

The  problem  with  the  speckle  noise  is  in  general  much 
smaller  by  direct  detection.  The  reason  for  this  is  that 
the  backscattered  field  has  a  finite  lateral  coherence 
length.  By  heterodyne  detection,  reception  takes  place 
over  only  one  coherence  area.  By  direct  detection,  on  the 
other  hand,  there  usually  is  a  spatial  averaging  over  many 
coherence  areas  in  the  receiver  aperture.  Thus,  the  vari¬ 
ance  of  the  intensity  fluctuations  is  correspondingly  de¬ 
creased  in  that  case.  Furthermore,  averaging  in  the  fre¬ 
quency  domain  also  may  take  place  in  the  direct -detec¬ 
tion  case. 

The  purpose  of  this  report  is  to  demonstrate,  in  more 
detail,  the  difference  between  direct  detection  and  hetero¬ 
dyne  detection  with  regard  to  the  speckle  noise.  This  is 
done  in  Chapters  2-4.  Then,  in  Ch.  6,  I  consider  the  sen¬ 
sitivity  limitations  imposed  by  the  speckle  noise  on  a 
range- resolving  differential  absorption  lidar  with  hetero¬ 
dyne  detection. 

2.  DIRECT  DETECTION 

I  shall  derive  an  expression  for  the  signal -to- noise  ratio 
which  takes  into  account  the  speckle  noise.  The  incident 
light  intensity  is  1.  The  detected  power  becomes 
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where  the  integration  is  over  the  receiving  aperture  area 
A.  The  mean  value  of  the  detected  power  is 

P0=</>>  =  Jj\/>d*r  =</>,!. 

There  are  two  noise  sources  by  the  detection:  noise  from 
the  detector,  and  noise  inherent  in  the  received  power  P. 
The  detector  noise  is  in  the  IE-spectral  region  usually 
thermal  noise  and  can  be  represented  by  a  noise-equiva¬ 
lent  power  NEP.  The  noise  inherent  in  the  received  power 
may  be  set  equal  to  aP,  the  standard  deviation  of  P.  Since 
these  two  noise  sources  are  independent,  the  total  noise 
becomes 

N=lf(NEFf+~i rj. 


Thus,  the  signal-to-noise  ratio  by  direct  detection  can  be 
written 


(WU- 


p» 

yiNtrf+o}' 


a) 


We  shall  now  calculate  the  variance  oJ=<P*>  — <P>*  in 
terms  of  the  statistical  properties  of  the  incident  light. 

where  /j-Ifo)  and  /*«!( r*). 

We  introduce  the  normalized  intensity  covarianoe  by 

0 ■  5 

and  assume  that  Cj  only  depends  on  £=>  |  rx  —  rt  ] ,  i.e., 
Ct^Cj(g).  Then, 

C,{g)#rt.  (2) 

Now,  the  backscattered  field  is  to  a  good  approximation 
a  complex  Gaussian  process,  i.e.f  the  amplitude  is  Ray¬ 
leigh  distributed,  the  phase  is  uniformly  distributed,  and 
the  intensity  has  a  negative  exponential  distribution.  This 
is  true  when  the  field  is  backscattered  from  an  incoherent 
target  (e.g.,  atmospheric  aerosol  particles)  and  in  the  ab- 
senoe  of  atmospheric  turbulence  (Goodman,  1976).  But  it 
is  also  true  both  at  sufficiently  weak  turbulence  and  strong 
turbulence  (Lee  et  al.,  1976;  Clifford  et  al.,  1978;  Pincus 
et  al.t  1978).  For  intermediate  values  of  the  turbulenoe, 
Pincus  et  al.  (1978)  experimentally  found  a  small  devia¬ 
tion  (at  most  about  12%)  from  the  unity  value  of  the 
standard  deviation  of  the  intensity  predicted  by  the  com¬ 


plex  Gaussian  process  model  of  the  backscattered  field. 
Thus,  in  the  following  we  assume  the  backscattered  field 
to  be  a  complex  Gaussian  process. 

Then  (see  Appendix  A), 


and 

where 


0*/ =  </>*. 

(3) 

CAe)-r*(eh 

(4) 

y(0)  »<*>**> 

m  <w> 

is  the  coherence  function  for  the  received  field  E.  As  be¬ 
fore,  EX*=E( r*),  etc. 

Substitution  of  Eqs.  (3)  and  (4)  into  Eq.  (2)  gives 


JMI 


It  is  convenient  to  define 


m 


*•  (5) 


Physically,  m  can  be  interpreted  as  the  number  of  correla¬ 
tion  cells  over  the  aperture  area  (Goodman,  1965a). 

We  obtain 


4= 


n 

m  * 


and  finally  from  Eq.  (1) 

(aiN)t--]=£j==.  («) 


The  expression  for  m  [Eq.  (5)]  can  be  simplified  in  the 
following  way:  Introduce  the  pupil  function 

(1  when  rt, 

\o  whenr>r#, 

where  r9  is  the  radius  of  the  circular  area  A.  Furthermore, 
change  the  variables  of  integration  by 

p-rj-r* 

Then 

m  -  (JJ  ^  y*(e)  JJ  ^  to)  D{*  -  Jp)  <***) 

where 
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Fig.  1 


It  is  not  difficult  to  realize  that  Jf(p)=the  shaded  area 
in  Fig.  1.  By  simple  geometrical  arguments,  we  find 

0  for  g  >  2rd, 

2ro  arcco8^--^^r5“  j  forp<2r0. 

Thus, 

r*(e)M(e)edej 

=  ( J*o  e (are  eos  e  —  g  ]/l  ~e*)yV*)<ie)  ,  (7) 

where  d  is  the  diameter  2rd. 

The  coherence  function  y(g)  appearing  in  the  expression 
for  m  has  recently  been  studied  by  Clifford  et  al.  (1978) 
and  Yura  (1978).  They  consider  a  Mar  system  with  back- 
scattering  from  aerosol  particles  and  take  into  account 
near- field  effects  and  atmospheric  turbulence.  The  field 
distribution  ET  over  the  lidar  transmitter  aperture  is  as¬ 
sumed  to  be  Gaussian  (this  corresponds  to  the  lowest  - 
order  mode  of  a  laser  with  confocal  mirrors): 

(8) 

Here,  E0  is  the  field  at  the  aperture  center,  r  is  the  dis¬ 
tance  from  the  center,  a  is  the  lie  intensity  radius,  k  = 
2njX,  and  /  is  the  focal  length.  The  coherence  function 
becomes 

y(e)»exp|-^)*J,  (9) 

with  the  coherence  length  qc  given  by 

(?)’  (s)’+  (o*k)  1  "•  W 

Here,  z  is  the  distance  to  the  scattering  region,  and  g0  is 
the  coherence  length  for  a  spherical  wave  which  has  prop¬ 
agated  a  distance  z  through  a  turbulent  atmosphere, 


Fig.  2.  (A)  The  number,  m,  of  correlation  cells  in  the  receiver 
aperture  as  a  function  of  rf/gc,  where  d  is  the  aperture  diameter, 
and  ge  is  the  lateral  coherence  length  of  the  incident  light.  (B)  The 
asymptote  i(d/pc)#. 

{>„=  "*• 

Cn  is  the  inde  x  -  of  -  refraction  structure  function  and  is  a 
measure  of  the  fluctuations  in  the  index -of-refraction. 
Substitution  of  Eq.  (9)  into  Eq.  (7)  gives 

^  ( j*o  e (»rc c°s e — e l^1 — e*) «*p [ ~ 2 (^)  ]  **e)  • 

(H) 

This  function  w  =  wi(d/pc )  was  calculated  numerically  and 
is  shown  in  Fig.  2.  Asymptotically  m  tends  to  4(d/pc)*, 
which  can  be  found  by  approximating 

p(arc  cos  p-p^l  -p*)  =  \nq 

in  the  integrand  in  Eq.  (11).  The  integral  can  then  be 
solved  analytically. 

To  demonstrate  the  effects  of  the  speckle  noise  on  a 
direct-detection  lidar,  we  now  consider  the  following  ex¬ 
ample  of  lidar  data: 

NEP=10~7  watt 

Transmitted  pulse  energy  =  I  J 

Optical  efficiency  =0.3 

Backscattering  coefficient  =  10~7  m_1  sr_l 

Atmospheric  attenuation  coefficient  =0.1  km-1 

Receiver  diameter  =0.5  m 

Radius  of  transmitted  beam  (the  quantity  a  in  Eq.  (8))  = 
0.1m 

/«oo  (collimated  beam) 

Wavelength  » 10  pm. 

Figs.  3-5  show  results  of  numerical  calculations  obtained 
from  these  assumptions.  Fig.  3  shows  the  coherence  length 
pc  as  a  function  of  the  distance  z  calculated  from  Eq.  (10). 
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Fig.  3.  (A)  The  lateral  coherence  length,  qc,  as  a  function  of  the  dis¬ 
tance,  *,  for  a  case  with  atmospheric  turbulence  (Cn  “  10-*14  m-8/4). 
(B)  Same  for  a  case  with  no  atmospheric  turbulence  (Cj  -  0). 

Results  are  shown  both  for  the  case  where  there  is  no 
turbulence  {Cl  =0)  along  the  propagation  path,  and  where 
there  is  turbulence  (<75  =  10”w 

Fig.  4  shows  the  corresponding  values  of  the  number 
m  of  correlation  cells  calculated  from  Eq.  (11). 

Finally,  Fig.  5  shows  (<S/JV)d  calculated  from  Eq.  (6). 
The  atmospheric  turbulence  does  not  affect  the  result,  at 
least  not  as  long  as  Cj  $  10~14  m“2/S.  In  the  figure  is  also 
shown  PJNEP ,  which  is  the  value  of  the  signal-to-noise 
ratio  without  regard  to  the  speckle  noise. 

It  is  seen  from  Fig.  5  that  there  is  a  considerable  de¬ 
crease  in  the  signal-to-noise  ratio  at  short  distances  from 
the  lidar  due  to  the  speckle  noise.  In  a  real  direct-detec¬ 
tion  lidar  system  this  decrease  can  be  smaller  for  a  number 
of  reasons. 

One  reason  is  that  in  addition  to  the  spatial  averaging 
just  discussed,  there  may  also  be  an  averaging  in  the 
frequency  domain.  The  necessary  wave-number  change 


Fig.  4.  (A)  The  number,  m,  of  correlation  cells  in  the  receiver 
aperture  as  a  function  of  the  distance,  *,  for  the  case  with  atmo¬ 
spheric  turbulence  (Ct  -  10*14  (B)  Same  for  the  case  with 

no  atmospheric  turbulence  (C&-0). 


Fig.  5.  The  signal-to-noise  ratio,  SjNt  as  a  function  of  the  distance, 
*.  (A)  Without  regard  to  the  speckle  noise  (PJNEP).  (B)  With 
regard  to  the  speckle  noise.  The  atmospheric  turbulence  does  not 
affect  the  results  of  this  example,  at  least  not  as  long  as  Cn  £ 


in  order  to  get  a  new  independent  sample  of  the  back- 
scattered  field  is  1//,  where  l=ct9  is  the  pulse  length,  c  = 
velocity  of  light,  and  tp  the  pulse  width  (Marshall  & 
Hitschfeld,  1953;  Wallace,  1953).  Assume  that  the  spec¬ 
tral  width  of  the  backscattered  light  is  Ak>\  jl  wave- 
numbers.  Then  the  number  mh  of  independent  samples 
over  which  averaging  occurs,  is 

mk  «  l\k.  (12) 

In  calculating  (S/N)d  from  Eq.  (6),  one  has  to  replace  m 
by  mmk. 

The  spectral  width  of  the  backscattered  light  may  either 
be  defined  by  the  initial  laser  line  width  or  by  the  Doppler 
broadening.  If  the  Doppler  broadening  dominates,  AJfc  = 
v/Ac,  and 


where  f,=A/v  is  the  time  scale  for  the  speckle  noise  due 
to  the  turbulent  velocity  v  in  the  scattering  volume. 

As  an  example,  for  a  NASA  direct-detection  DIAL 
under  development  (Stewart  &  Bufton,  1978),  the  laser 
line  width  will  be  0.033  cm-1  and  the  pulse  width  50  ns. 
Then,  from  Eq.  (12),  m*=50. 

The  speckle  noise  can  of  course  also  be  reduced  by  aver¬ 
aging  in  the  time  domain.  If  the  integration  time  or  gate 
time  <fl>Min(fp,  t9),  then  there  will  be  an  averaging  over 
mt  independent  samples,  where 

t9 

Tn,'=  Min 

A  further  reason  for  the  signal-to-noise  decrease  to  be 
smaller  than  shown  in  Fig.  5  is  that  the  spatial  averaging 
over  the  receiver  aperture  can  be  more  effective  than  cal¬ 
culated.  In  the  derivation  of  Eq.  (10)  it  was  assumed  that 
the  laser  was  in  single-mode  operation.  However,  for  a 
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laser  in  multi-mode  operation,  the  beam  divergence  can 
be  larger.  This  results  in  a  smaller  coherence  length  gct 
a  larger  number  of  correlation  cells  over  the  receiver  aper¬ 
ture,  and  a  more  effective  spatial  averaging.  If  atmo- 
spheric  turbulence  effects  are  negligible,  the  coherence 
length  can  be  estimated  from  the  van  Cittert-Zemike 
theorem  (Boro  &  Wolf,  1965,  p.  508).  Approximately,  in 
the  far  field  the  coherence  length  qc  becomes  A/0,  where  0 
is  the  beam  divergence. 

Finally,  1  want  to  point  out,  that  in  this  analysis  at¬ 
mospheric  turbulence  has  been  considered  only  with  re¬ 
gard  to  its  effects  on  the  coherence  length  ge.  But  the 
turbulence  also  gives  rise  to  extra  noise  (scintillations)  in 
the  signal.  This  effect  has  not  been  taken  into  account. 
However,  it  is  probably  often  a  small  effect,  as  mentioned 
before. 


3.  HETERODYNE  DETECTION 


In  the  same  way  as  Eq.  (1 )  was  used  for  a  direct-detection 
system,  we  set  for  the  signal-to-noise  ratio  of  a  lidar  with 
heterodyne  detection 


(SlN)k^ 


Po 

l/iNEPf- t7*#’ 


(13) 


P0  is  now  no  longer  the  intensity  integrated  over  the  total 
receiver  area,  but  instead  (Goodman,  19656) 


with  the  effective  area 


A‘-i!S/,~SLr,er'tA- 


There  is  no  spatial  averaging  of  the  intensity  fluctuations. 
Therefore,  assuming,  as  before,  a  negative  exponential 
distribution  for  the  intensity, 

aP  =  P0. 


Op  will  be  decreased  by  the  factor  J3*p.  From  Eq.  (13) 
we  get 


It  is  easy  to  show  that  the  highest,  in  this  way,  attainable 
signal-to-noise  ratio  is 

(S/tf)*“0.89[(S/tf  ).]*'*, 

where 

is  the  signal-to-noise  ratio  in  the  absence  of  speckle  noise 
and  with  2?—  l/f„. 

The  corresponding  optimal  value  of  the  IF  bandwidth 

fl_0.79[  (g/W.  (M) 

fp 

A  necessary  condition  for  these  results  is  of  course  that 
the  bandwidth  of  the  backscattered  light  is  ^  B. 

4.  COMPARISON  BETWEEN  DIRECT  DETECTION 
AND  HETERODYNE  DETECTION 

It  may  be  of  interest  to  try  to  compare  a  lidar  system 
with  direct  detection  and  one  with  heterodyne  detection. 
For  simplicity,  let  us  assume  that  the  two  systems  are 
identical  (same  pulse  energy,  optical  efficiency,  etc.)  ex¬ 
cept  for  the  detection  technique.  The  noise-equivalent 
power  for  the  direct-detection  system  is  (NEP)d  and  is 
determined  by  the  detector  thermal  noise.  For  the  hetero¬ 
dyne  case  it  is  ( NEP)h>  the  shot  noise  induced  by  the 
local  oscillator.  Typically,  at  IR  wavelengths  (NEP)d  is 
larger  than  (NEP)h  by  several  orders  of  magnitude.  How¬ 
ever,  in  the  direct-detection  system  the  speckle  noise  is 
effectively  averaged  out  in  the  spatial  and  frequency  do¬ 
mains,  and  it  may  often  be  possible  to  obtain 


The  NEP  now  is  the  noise-equivalent  power  due  to  the 
local  oscillator  shot  noise, 

NEP-—. 

V 

Here  h  is  Planck’s  constant,  v  the  frequency,  B  the  IF 
bandwidth,  and  r)  the  quantum  efficiency.  Often,  P0> 
NEP,  and  then  (8/N)h  ~1. 

Usually  one  chooses  Since  the  necessary  fre¬ 

quency  change  to  get  an  independent  sample  of  the  back- 
scattered  field  also  is  l/tp,  we  see  that  there  can  usually 
be  no  averaging  of  the  speckle  noise  in  the  frequency 
domain. 

However,  if  P9>NEPt  we  can  in  fact  increase  the 
signal-to-noise  ratio  by  choosing  £>!/<*.  The  variance 


In  the  heterodyne  case,  on  the  other  hand,  we  get,  as 
was  shown  above,  at  most 


Therefore,  the  direct-detection  technique  gives  a  higher 
signal-to-noise  ratio,  at  least  as  long  as 


or 


p j 

(NEP), 


>0.89 


P,  >  0.86 


(NEPft* 

(nep) r 
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If  the  IF  bandwidth  is  not  optimized  in  the  sense  of 
Eq.  (14),  but  is  B  =  l/fp,  we  get  instead 

P,>(NEP)+ 

In  conclusion,  the  direct-detection  system  has  the  largest 
aignal-to-noise  ratio  at  short  distances,  whereas  at  longer 
distances  the  heterodyne-detection  system  is  superior. 


5.  DIFFERENTIAL  ABSORPTION 

In  this  chapter  we  will  consider  the  sensitivity  limitations 
imposed  by  the  speckle  noise  on  differential -absorption 
measurements  with  a  heterodyne-detection  lidar. 

To  start  with,  let  us  set  up  the  basic  equations.  The 
lidar  equation  can  be  written 

p»  =  exP  [  -  2  Jo  +  a)  (fcj ,  (15) 

where  P0  is  the  average  received  power,  /9(z)  is  the  back- 
scattering  coefficient  at  distance  z,  cr  is  the  absorption 
cross  section  for  the  gas  of  interest,  a  is  the  attenuation 
coefficient  due  to  all  other  constituents  of  the  atmosphere, 
Ng  is  the  gas  concentration  at  distance  z,  and  C  is  a  con¬ 
stant. 

It  follows  from  Eq.  (15)  that  the  average  gas  concen¬ 
tration  N0  over  the  distance  z1  —  z2f  i.e., 


is  given  by 


where 

r  =2[ff(A1)-<r(A8)](z1-22), 

o^)  =  molecular  absorption  cross  section  on  an  absorp¬ 
tion  peak  of  the  gas  of  interest, 

=  molecular  absorption  cross  section  off  the  absorp¬ 
tion  peak,  and 

P0(A,,  zf)  =  averaged  received  power  at  wavelength  A*  and 
range  z}. 


j: 


N2dz 


-H 


1  1  P p(^2»  2l)  Zz) 

0  r  PAKzMKztY 


(10) 


For  the  derivation  of  Eq.  (16)  it  was  assumed  that  Xx  and 
Aa  are  so  close  to  each  other  that  the  backscattering  co¬ 
efficient  ($  and  the  attenuation  coefficient  a  have  the  same 
values  at  the  two  wavelengths.  This  is  a  quite  reasonable 
assumption  for  fi,  but  may  not  always  be  true  for  a  be¬ 
cause  of  interference  from  other  gases.  However,  by  dif¬ 
ferent  compensation  techniques  these  interference  effects 
can  be  made  negligible  (Murray,  1978). 

An  estimate  N  of  the  gas  concentration  is  obtained  by 


V  M  1  In  *  *t) 

T  P(XX9Zl)P(XtyZt)' 


(17) 


We  will  consider  the  large-signal  case,  i.e.,  where  Pe> 
NEP.  Then  the  uncertainty  aN  in  the  estimate  N  is  due 
to  fluctuations  in  the  received  power  P.  These  fluctua¬ 
tions  in  turn  are  due  to  speckle  noise,  atmospheric  tur¬ 
bulence,  and  variations  in  the  backscattering  coefficient 
p  and  the  attenuation  coefficient  a. 

The  time  scale  for  the  turbulenoe-induoed  fluctuations 
is  I^Az/t?,  where  v  is  the  wind  velocity  transverse  to  the 
propagation  path  (Clifford  et  ah,  1978).  Typically,  Vkzjv  ** 
10~*  — 10_*  s.  The  time  scale  for  the  variations  in  back- 
scattering  coefficient  ft  and  the  attenuation  coefficient  a 
is  larger  than  10~3s  (Schotland,  1974).  It  is  therefore 
strongly  to  recommend  that  the  measurements  on  and  off 
the  absorption  line  be  made  within  a  time  which  is  smaller 
than  about  10-3  s,  since  then  the  fluctuations  due  to  tur¬ 
bulence  and  variable  aerosol  structure  cancel  out  in  the 
power  ratios  in  Eq.  (17).  The  best  is,  of  course,  if  the 
power  measurements  are  made  simultaneously  (Derr 
et  ah,  1974;  Stewart  &  Button,  1978). 

The  wavelength  dependence  of  the  turbulence  effects 
in  Eq.  (17)  is  negligible.  The  characteristic  transverse 
length  for  correlation  of  the  intensity  fluctuations  is  VJz 
(Fante,  1975).  Furthermore,  in  a  differential-absorption 
measurement  the  fractional  wavelength  shift  typically  is 
about  0.2%.  Therefore,  the  power  scintillations  at  Ar  and 
A2  can  be  considered  to  be  completely  correlated,  and  they 
cancel  out  in  the  power  ratios  in  Eq.  (17). 

In  the  following  we  assume  that  the  measurements  on 
and  off  the  absorption  line  are  made  effectively  simul¬ 
taneously.  The  only  contribution  to  the  uncertainty  aN 
then  comes  from  the  speckle  noise.  The  four  samples 
,P(A2,  zx),  P(Ax,  z2),  P(AX,  zL),  and  P(A2,  z2)  are  certainly  in¬ 
dependent  with  regard  to  the  speckle  noise.  This  follows 
since  an  independent  sample  of  the  speckle  noise  is  ob¬ 
tained  if  \zx  —  z2|  >1  or  | Aa  —  A2)  >A2/Z,  where  l^ctp  is  the 
pulse  length. 

From  Eq.  (17)  we  get 

Var(2V)  =  ^  {Var[In  P(Xi,  z,)]  +  Var[ln  P(A„  z,)] 

+  Var  [In  P(A„  2l)]  +  Var[ln  P(^,  *,)]}, 

where  Var(-)  denotes  the  variance.  The  probability  dis¬ 
tribution  function  for  P  is 

WH-±«p(-£ ), 

where  as  before,  denotes  the  mean  value.  We  obtain 
(Gradshteyn  A  Ryzhik,  1965,  4.331  and  4.335): 

Var  (In  P)  -  (ln*P>  -  <ln  P># 

-  J*p(P)  In *PdP-  p(P)  In  PdPj* 
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Fig.  6.  The  standard  deviation,  ton ,  of  the  estimate  of  the  gas 
concentration  as  a  function  of  the  number,  m,  of  independent 
samples.  The  function  is  calculated  from  (A)  an  approximate  ex¬ 
pression  (Eq.  (19)),  and  from  (B)  an  exact  expression  (Eq.  20)). 


Thus, 

2rr2 

Var(V)=— 2) 

and 

i/it — 2.6 

%=KVar  (#)«_«—. 

ri/6  r 

In  a  sensitivity  analysis  it  is  reasonable  to  set  the  minimum 
detectable  gas  concentration,  Nmln ,  equal  to  the  standard 
deviation  oN.  Thus,  Nf03n  =  2.6fT.  If  we  average  over  n 
pulses,  the  standard  deviation  of  course  instead  becomes 
OxIVn. 

Finally,  it  can  be  of  some  interest  to  consider  the  case 
where  there  is  averaging  of  the  speckle  noise  in  the  fre¬ 
quency  domain,  as  discussed  in  Ch.  3.  Then  we  have  in¬ 
stead  of  Eq.  (17): 

j  2’  2>)  ,jqv 

with 

Q$i,  zj)  —  —  2  Pfii*  zj)> 

7tl  I 

and  m  is  the  number  of  independent  samples  of  P  over 
which  averaging  takes  place.  It  follows  that 

Var  (jV)  =  ^  {Var  [In  Q(^  2,)]  +  Var[ln  Q(A„  2,)] 

+  Var  [In  Q(AX,  2,)]  +  Var  [In  Q(A„  2,)]}. 

Now,  if  m  is  sufficiently  large,  then  the  probability  dis¬ 
tribution  for  Q  is  sufficiently  narrow,  and  we  can  approxi¬ 
mate  (Bevington,  1969,  p.  69): 


where  Q0*P0  is  the  mean  value  of  Q.  But 
Var  (P)  P% 


Thus, 


and 


Var(G)* 


m 


m 


Var  (In  Q)~~t 
m 


Var^^, 


T/m 


(19) 


It  is  possible  to  show  (see  Appendix  B)  that  for  a  general 
value  of  m  >  1 ,  we  instead  get 


2  )/n2  V  1 


(20) 


The  standard  deviation  calculated  from  Eqs.  (19)  and  (20) 
is  shown  as  function  of  m  in  Fig.  6. 
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APPENDIX  A 

RELATION  BETWEEN  INTENSITY  COVARIANCE  AND 
COHERENCE  FUNCTION 


The  field  E  is  assumed  to  be  a  complex  Gaussian  process, 
i.e., 

E  =*  iy, 

where  x=\E\  cos <j>  and  y=\E\  sin <f>  are  independent, 
zero  mean  Gaussian  random  variables  with  variances  a2. 
The  probability  distribution  is 


The  intensity  is 
Thus 


(^|2  =  x2  +  y2. 


<J>  =  <*2>+<l/2>=2a2, 


</*>  -  <**> + <y  >  4-  2<*2>  <y2>  ==  804. 


It  follows  that 

0?  =  </*>-</>*  =  </>*. 

vhich  is  Eq.  (3). 

We  shall  now  relate 


c  </,/,)-</>» 

J 


to  y 


<W>  ' 


</,  /,)  -  <(*f  +  y?)  (*S  + |4)> 

=  <*?*!>  +  <*f>  +  <*»>  <y?>  +  <3&yi> 

-2«*J4>+A 


where  we  have  used  that  <«*>  =  <t/2>  ^<ra  and,  by  sym¬ 
metry,  <y?yf>  *=<rfa|>.  {rfri)  can  be  expressed  in  <rjX2) 
in  the  following  way  (Burdic,  1968,  p.  277): 

Introduce  a  new  variable  z  by  x2  —  cx2  4-  z,  where  the 
constant  c  is  defined  by  c  =  (x1x2y/o2.  It  is  then  easy  to 
show  that  xx  and  z  are  uncorrelated,  i.e.,  (x1z>~0.  But 
since  xx  and  z  also  are  Gaussian,  it  follows  that  they  also 
are  independent.  We  then  find  that 

<*l4>  =  +  *)*>  =  3  <*!*,>*  +ff*<22). 

Furthermore, 

<**>-«*,—<*,)*>-**— 

Therefore, 

<44>  =  2  <*,*,>* +<x*. 

Now, 

<i>*  ** 

7  <|J?|*>  a*  • 

Combination  of  the  last  two  equations  finally  gives 

<w*, 

which  is  Eq.  (4). 
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APPENDIX  B 


DERIVATION  OF  A  FORMULA  FOR  THE  STANDARD  DEVIATION 
OF  ESTIMATED  GAS  CONCENTRATION 


The  basic  equation  is  Eq.  (18): 


with 


®  jjjjs*  h)Q$v  *i) 


Q$t,  */)  -  m  2  P(^i>  */)• 


(B-l) 


In  this  way  pN(N)  can  be  calculated.  The  result  is 
pw(y)°  [[(rT^l)!]8]  TexP[wT(Ar-A,“)] 


i 


t^di 


The  probability  distribution  pP(P)  for  the  received  power 
P  is  the  negative  exponential, 


From  this,  we  can  derive  the  probability  distribution 
pN(N)  and  the  standard  deviation  aN  for  the  estimate  N 
defined  in  Eq.  (B-l). 

First  we  note  that  the  probability  distribution  function 
Pq(Q)  for  Q  is  the  gamma  distribution  (Marshall  &  Hitsch- 
feld,  1953;  Wallace,  1953): 

where  Qq~P0  is  the  mean  value  of  Q  and  P. 

Next,  introduce  random  variables  Xl*  and  by 

Y  —  z|)  y  —  9^!’.??)  y  _  Y  Y 

1  Q{h,*iY  2  Q(KhY  3  125 


so  that 


r  * 


o  (t+lfm{t  +  exp[T(N-Xom^‘ 

Here  N  =  N0  is  the  mean  value  of  Nc,  and  is  defined  in 
Eq.  (16).  The  fcth  moment  of  N  becomes 

<N*>  =  (°°Nkp^N)dN 
*  J  -00 

\(2m~l)\y  f®  t2m~l  . 

=  [((m-l)!fj  rJ0  (t+l)2m<* 

•+00  ^6Xp[mT(JV-^0)] 

a/V  . 


/: 


!  -00  {*  +  exp  [r(X  -  X0)]}2m 

Change  of  variable  of  integration  in  the  last  integral  by 

exp[rW-N0)] 

t 

gives 

f  (2m—  1)!  I2  1  r  l  , 

<  *>  =  Lt(»-  1)!]2J  t*Jo  ( 

r  {In  [st  exp  (Tff„)]}*sM-1 

Jo  («+i  r 

If  we  furthermore  utilize  that 

{In  [st  exp  (rN0)]}*  —  (In  s  +  In  t  +  xN^)k 


Now,  in  general,  if  x  >0  and  y>  0  are  independent  random 
variables  with  probability  density  functions  px(x)  and 
pv(y),  and  if  z-xjyy  then  the  probability  density  function 
for  z  is 

1*00 

?>.(*) -I  yp*(y*)p»(y)<*y- 

If,  instead,  z  =  xyy  then 

P.(*) =  y  P*  (y)  P»(y) 

With  these  formulas,  the  probability  density  functions 
for  Xv  Xr  and  Xs  can  be  derived.  Furthermore,  the 
probability  density  functions  ps(N)  and  Px,(X8)  for  N 
and  X9  are  related  by 
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then 


<",>-?[t^iw]'.4  k  C)  0 ^  -w ' 


where 


Jo  (<+l)2ra  *' 


(B-2) 


The  problem  of  evaluating  <N>  is  thus  reduced  to  evalu¬ 
ating  integrals  of  the  type  in  Eq.  (B-2).  It  can  be 
shown  that 

Ap  *=  0  when  p  is  odd, 

[Jm~  1)1J» 

0  (2m- 1)1’ 

.  c(m-i))]V^  yn 

*  (2m- 1)!  \«  ,f,w 


i 


With  these  results  it  follows  that 


But  (Gradshteyn  &  Ryzhik,  1965,  8,360  and  8.363) 


<N>-N0f 


1 


tZo  (?+%)*  dz 2 


In  I», 


and 


where  T(^)  is  the  gamma  function. 

Also, 

In  T(z)  =  x  In  z  when  x-*>  oo. 


which  is  the  desired  result.  For  this  to  be  consistent  with  Therefore, 
Eq.  (19),  it  has  to  be  shown  that  aN  given  by  Eq.  (B-3) 
tends  to  2 /(r/m)  when  oo.  To  this  end,  we  first  note 
that 

®  1  7T* 

rt*!  P2"  6  1 

and  therefore 


1 


1 


,r0(v  +  *)2  x 


when  z~+  oo, 


and 


a  m-l  1  0C  1  oo 

^  X  ^  X 


1 


6  ,fir2  ,rQ(v  +  mf 


7*  "-1  1  1 

—  —  >  -4“*—  when  m-*- 00, 

6  ,fi  va  w 


—  — r=  when  m  -►  00 . 
Tr  w 


10 


FOA  Reports,  Vol.  13,  No.  1,  1979 


